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Abstract 



For top quark decays into heavy quark mesons T and B* , a com- 
plete calculation to the leading order both in QCD coupling con- 
stant ot s and in v, the typical velocity of the heavy quarks inside 
the mesons, is performed. Relatons between the top quark mass 
and the decay branching ratios are studied. Comparion with the 
results which are obtained by using the quark frangmentation func- 
tions is also discussed. The branching ratios are consistent (within 
a factor of 2 ~ 3) with that obtained using fragmentation functions 
at m t ~ 150 GeV. 
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I Introduction 



The success of Standard Model (SM)|T|]0|| suggests that the top quark must exit ||]. 
Recently, from the direct search at the Tevatron, the CDF and DO groups confirmed the 
existence of heavy top quark ||, which has a mass of (176 ± 8 ± 10) GeV or (199l 2 i ± 22) 
GeV. Then next the experimental studies will focus on the determination of its properties. 
Among them, the precision measurement of the top quark mass and of the production cross 
sections and distributions will certainly be in the first studies of interest. Since the fermion 
mass generation can be closely related to the electroweak symmetry breaking, one expects to 
find some residual effects of the breaking in accordance with the generated mass especially 
in the much heavier top quark sector. Therefore, it is important to study the top quark 
system as a direct tool to probe new physics effects^]. 

As the new found top quark mass is much heavier than that one previously expected, its 
decay widths in many processes should be reevaluated. With future upgrade and increases 
in total integrated luminosity, experiments at the Fermilab Tevatron will carry a carefull 
investigation of the production and decay of top quark. With the operation of CERN Large 
Hadron Collider (LHC), one expects to obtain roughly 10 7 — 10 8 ti pairs per year ]7J and 
then one may be able to observe various top decay channels, which will give further tests of 
the SM and heavy meson production mechamisms. On the other hand, the study of heavy 
quarkonium and heavy meson production in turn provides a ground to precisely test the 
perturbative quantum chromodynamics (PQCD). 

Within past few years, much progress has been made in the study of production mech- 
anism of heavy quarkonium and heavy mesons with large transverse momentum at high 
energies. The parton fragmentation in the heavy quarkonium and heavy meson production 
has been carefully studied || 1| and the universal fragmentation functions for different spin 
and angular momentum states have been obtained. With these functions, the branching 
ratios of heavy quarkonium and heavy meson production for various channels may be easily 
obtained if the fragmentation mechanism is in dominance. Based on the new factoriza- 



tion formulas[10] and the production mechanism, one may explain some new experiment 



results|ll]. Hence, in what energy region and in what degree of precision the fragmentation 
mechanism works in a specific process should be concerned. 

The most important concequence of a heavy top quark is that, because its lifetime is 
short, after created in free state it does not have time to bind with light quarks to form a 
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bound state[13|. Therefore, the final frontier for the study of heavy- quark- ant iquark bound 
states may be the B c and B c , which carry flavors explicitly and are the ground states of 
bound bc and cb systems. The properties of these mesons and the possibility to find them 
at existing accelerators have been discussed in Refs. H (14| . The top rare decay to B* is also 
an important channel to study the still undiscovered heavy quark-antiquark bound state at 
LHC in the near future, though it is not practical at present with limited top quark events. 

In this paper we present the full leading order calculation of top quark decays into B* 
and T in strong coupling constant a s as well as in v 2 , where v is the typical relative velocity 
of the heavy quarks inside the bound state. The reliance of the branching ratios of heavy 
quarkonium and heavy meson production in top decay on the mass of top quark is also 
studied. 



II Formalism 



The amplitudes for t — > W + cB* and t — > W + bT involve two Feynman diagrams respec- 
tively as shown in Fig.l. They can be writen down using standard Feynman rules for bound 
state productions and decays [12|] |19| or from the description of Bethe-Salpeter equation for 
bound states. In the following we will start with discussion from the latter concept. 



B* C (T) 




V B*(T) 




Fig.l 



Based on the Mandelstam formalism [16] the amplitude of top decays corresponding to 
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diagrams Fig.l is written as 

M = ^Lu(q 2 ,s) J ^{A 1 + A 2 }v( qi ,s). (1) 

Here s, s are the spin projections of the quark and the antiquark; qi and q 2 are their momenta; 
g s is the coupling constant of QCD; g is the weak coupling constant in WS model. 
The general form of the partial width for the top decays under consideraton is 

r« threeiody) = £ <fo 2 ds^\M\ (2) 

In the following we will discuss the two channels of top quark decaying into vector mesons 
B* and T which might be in the first place being measuerd in the future experiment on top 
rare decay. 

A. B* + W + + c 
At leading order in a s , the corresponding Ai,A 2 in Eqn.l for t — > B* + W + + c are 

Al = ^771 \2'~m \2 2 ?w P L, (3) 

W + P2) \P + Vi) - 

A 2 = Jf M> < ^"^V. (4) 
\P + P2) \Px + «) - mi 

Where the Pl = |(1 — 75); x{l) is the Bethe-Salpeter(BS) wave function of the5* meson 

with inner relative momentum q between the heavy quarks; t w is the polarization vector 

of W + boson. In Eqn.(|3[)(|]), p,k,p',pi and p 2 are the momenta of B* meson, W + boson, 

c quark, b quark, and c quark, respectively. For B* meson, the two quarks have different 

masses inside the bound state. Their momenta p\ and p 2 have the relation 



Pi = ViP + <? > P2 = V2P - q, 



(5) 



where rji 



and rj 2 



-. Under the instantaneous approximation [H| with the 



(m h +m c ) 1 1 (m b +m c ) ' 

negative energy projectors being neglected, the BS wave function of the bound state 
may be expressed as 



x{q) 



M — E 1 — E 2 



2tt {P 10 -E 1 )(P 20 -E 2 



Mo)- 



(6) 
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Here M is the mass of bound state; Pio and P20 are the time components of quark and 
antiquark momenta inside the meson, and Ei,E 2 are their energies. From the standard BS 
wave functions in the approximation that the negative energy projectors are omitted, the 
vector meson wave function can be projected out as : 

$ (<?1 = ^E< JM l lsL ™> A U<?>0 AM+ P) lQ M{-q)^ Lm {q), (7) 

sm 

where e is the polarization vector associated with the spin triplet states. h} + (q) and h 2 _(—q) 
are positive energy projection operators of both quark and antiquark . 

m = E+1 °\l +mc '<\ w-$= E+ '**2*- nv » . (8) 

After taking the nonrelativistic approximation the bound state wave function may be 
reduced to a simple form. As in Refs. |i8|jl9|| , we also ignore the dependence of A\ and A 2 



on the relative momentum q. It may be considered as the lowest order approximation. After 
considering the normalization of $(<f), we can get 



Here R§*(0) is the radiual wave function at origin. Hence the matrix element squared can 
be carried out directly, 

^|M| 2 = Y,(\M 1 \ 2 + \M 2 \ 2 + 2ReM 1 M;). (10) 

The square of the full amplitude is complicated and lengthy, it can be found in Appendix A. 
The kinematical variables entering Eq.(Q) are defined as 

Sl = (p'+p) 2 , S2 = ( p+A; ) 2 , (11) 

in which the phase space boundaries are readily found to be 



and 



mi+M^- — [(s 2 -mi+mi)(s 2 +M 2 § ,-mi)TXHs2,mlml)X^s 2 ,M^,ml)], (12) 

AS 2 



s 2 = (Ms. + m w ) 2 , (13) 
s+ = (m t - m c ) 2 , (14) 
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where 



\(x, y,z) = (x - y - zf - Ayz. 



(15) 



B. t^T + W+ + b 

The A\ and ^4. 2 in Eq.(|TJ) for t — > W + bT are the same form as Eqs.(|3])(f|) except x{l) 
representing bottomonium BS wave function and c(c) replaced by b(b) quark. Under the 
same argument as in A, we have 



r d^a 1 

I^Wf^4Vm AF + Mr)Rm - <16) 



The matrix squared is given in Appendix B. The phase-space bounderies to be used in the 
numerical evaluation of Eq.(|2|) in this channel are given by 



sf = ml + m\- ^—[(s 2 -m 2 t +ml)(s 2 + m 2 r -ml)^ ^H s 2,m 2 t ,m 2 b )X^(s 2 ,rn\,m 2 w )], (17) 
zs 2 

and 



s 2 = (m T + m w ) 2 , (18) 
4 = (m t -m b ) 2 . (19) 

III Numerical Calculation and Results 

We can now make use of this formalism to evaluate the decay rates of top quark to T 
and B* mesons. In numerical calculations, we take the parameters as follows: 
m b = A.9GeV, m c = 1.5GeV, m t = 176GeV, m w = 80.22GeV. 
a s = 0.26 for t -> W + cB* c , 

a s = 0.19 fort-*w + bT. (20) 

Under the non-relativistic approximation, 

M T = 2m b , M B * = m b + m c . (21) 

The radial wave functions at the origin for B* and T can be estimated through potential 
models]!?]] and electric decay rate. As in Ref.|| we take 
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R s .(0) 2 = (1.18GeVf, 

\Rr(0)\ 2 = (l.8GeVy j . (22) 
Finally, we obtain the partial width of the decay t — > B* + W + + c : 

T(t -> B*W + c) = OMMeV, (23) 

giving 



T(t -> 5W + c , 
; = 4.12 x 1(T 4 . 

r(t -> w+b) 

The partial width of the decay t —>■ T + iy + + 6 is 



(24) 



r(t -> T + W + + 6) = 1.54 x 10- 2 MeV, (25) 

giving 

T(F^W = °' 98 x 10 (26) 

The complete leading order calculation of the decay rate for t — > W + bT has already 
been calculated [0] for top quark with a mass of 100 GeV and slightly different values for 
other parameters, giving the branching fraction 4 x 10~ 7 . In our calculation , we find 
the branching ratio is two order of magnitude greater when the top quark mass is set to 
176 GeV. This indicates that the T and B* meson productions in top decay is really a 
quark fragmentation process if the top quark is heavy enough. Using the fragmentation 
function obtained in Ref.||, one can readily obtain the branching ratios for t — > B*W + c and 
t — > TW + b. However, in how much degree the univeresal fragmentation functions suit the 
top rare decays is not yet very clear. To answer this question, we give out the relationships 
between the top quark mass and the branching ratios in Fig. 2 and Fig. 3 respectively for 
t — ► TW + b and t — > B*W + c. From Fig. 2 and Fig.3 we can see that the branching ratios 
increase rapidly as the top quark mass increases when m t < 120 GeV, and the heavier the 
top quark is, the less the branching ratios change. When top quark mass is heavier than 
200 GeV, the branching ratios seem to be saturated and no longer depend on the top quark 
mass. That means the univeral fragmentation functions works well. As for the numerical 
results, with the same parameters as used in Refs.0 at m t = 176 GeV, our complete leading 
order calculations give R(t — > TW + b) = 0.98 x 10~ 5 , which is about two times smaller than 
in Ref.|J using the universal fragmentation function and R(t — > B*W + c) = 4.12 x 10~ 4 , 
which is almost as large as Ref.|| using the universal fragmentation function. 

In Fig. 4 and Fig. 5 we also show the relations of these branching ratios to the strong 
coupling constant a s . 
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IV Conclusion 



We have presented two dominant three-body decays of top quark to ground state heavy 
mesons. The decay widths and branching ratios are O.QAMeV and 4.12 x 1CT 4 for B* 
productionin ; 1.54 x 10~ 2 MeV and 0.98 x 10" 5 for T production . As about 10 8 it pairs per 
year will be produced at the planned LHC, this means it is possible to accumulate about 10 4 
B* events and 10 3 T events a year. This will make it possible to discover the B* meson and 
to further study its properties experimentally. At NLC (Next Linear Collider), rare decays 
of top quark could be searched for down to a very small branching ratos, and a systematic 
and general search for top rare decay modes may be possible. Then the two top decay 
modes discussed in this paper might correspond to a detectable level at NLC. On the other 
hand the results obtained in this paper will also be a helpful reference to other production 
mechamisms concerning the top quark decay. 
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APPENDIX A 



Following we give the expressions of squared matrix for t — > W + cB*. The vertex co- 
efficents are not taken into consideration. In this appendix we take M§* = §Mg» just for 
calculation convenience. 



{Am\M%,m 2 w + %m b M% t m c m 2 w + AM% t m 2 c m 2 w ) 

JD C D c JD C 

1 1 

(771S2 - 4tt2M|. + r) 2 ml + 4r/ 2 M|, - m 2 ) 2 4^ - m 2 ) 2 

x {<6Am 2 b M%,(k ■ p) 2 (p ■ p')(t ■ p) - Mm b M 2 B ,m 2 t (m b + m c )(k ■ p) 2 (p ■ p') 

+ 128m 2 M| ; m c (£; • p) 2 {t ■ p) - 320m b M^m c m 2 (m b + m c )(k ■ p) 2 

+ 32(4m 2 M|* + m 2 m 2 + 4m 2 m 2 + 2m b m c m 2 

+ 8m(,m c m^ + m 2 m 2 + 4m 2 m 2 u )M|„ (k ■ p) 2 (t ■ p') 

— 64Mg t (2m 2 Mg* — rn\m 2 — m 2 h m 2 w — 2m b m c m 2 

— 2m b m c m 2 w - m 2 m 2 t - m 2 m 2 w )(k ■ p)(k ■ p') (t ■ p) 

+ \2%m\M% t {k -p)(k- t)(p ■ p) + %Am b M% t m 2 w (m b + m c )(k ■ p)(p ■ p')(t ■ p) 

+ 64M|» (— 4m 2 m c M|» + m 2 m 2 m c + Am 2 w m 2 m c + 2m 2 t m 2 c m b 

+ %m 2 w m b m 2 c + m^m 2 + Am\m 2 w ){k ■ p)(k ■ t) 

— lQ0Ml i m 2 m 2 w (m b + m c ) 2 (k-p)(p-p') 
+ l§2m b m c Mg«m 2 w (m b + m c )(k-p)(t-p) 
+ 704m b M^,m 2 w (m b + m c )(k ■ p)(t ■ p') 

— 672m c Mg*m 2 m 2 v (m c + m b ) 2 (k ■ p) 

— 64M|» (4m 2 M|„ + m\m 2 t — 2m\m 2 w + 2m b m c m 2 — Am b m c m 2 

+ m 2 c m 2 — 2m 2 m 2 w )(k ■ p')(k ■ t) + l92m b Mg*m 2 u (m b + m c )(k ■ p')(t ■ p) 

— 128M| *m 2 t m 2 w (m b + m c ) 2 (k ■ p) + 7§&m b m c M\,m 2 w (m b + m c )(k ■ t) 
+ l92m b M^m 2 w (m b + m c )(k-t)(p-p') 

+ 16M|,m 2 (8M| ; m 2 + 2m 2 m 2 - m^m^ 

+ Am b m c m\ — 2m b m c m 2 u + 2m 2 c m\ — m 2 c m 2 w )(p ■ p')(t ■ p) 

— h\2m b M\,m 2 m 2 w (m b + m c )(p ■ p) 

+ 64M| ; m c m 2 (8m 2 M| ; + 2m 2 m 2 - m^m^ 
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+ 4m b m c m 2 — 2m b m c m 2 w + 2m 2 m 2 — m 2 m 2 )(£ • p) 

+ 32M|.m^(28mgM|. + 7m 2 b m 2 t 

— hm^rr? w + YIm b m c m\ — \§m b m c vr? w + 7m 2 c m\ — 5m 2 m 2 )(t • p') 

— 1408m b m c M^rn 2 m 2 u (m b + m c )} 



4?7l(si - m^) 2 (si - m 2 ) 2 
x {(8M§, -m2 + mg)(p-p')(*-p) 

+ 2M|» (-4M|» - 8M s *m b + 8M s *m c + m 2 c + m 2 b - <6m b m c )(t ■ p') 
+ 4M|, (m 2 - 3m 6 m c )(t • p) + 8M§. (p • p')(* * p') + 2 (p • p') 2 (* • V ) 

+ =-£-(-4M|. + m 2 >m 2 ^6m f) m c )(/c-t)(/c-p / ) 

16M? 9 ~ 
+ ^ (A; • -p^ip-p') H 5- (8Ml ; + m 2 - m 2 c ) ( fc • t) (k ■ p) (p • p') 



8M? 4 
+ — (m 2 - 3m fe m c ) (k ■ t) ( fc • p) + — ( k ■ t) ( fc • p') (p • p') 2 ) } 



1 1 1 

M m m 2 w (m b + m c ) («i - ™ 2 ) 4?$ (si - m 2 ) 2 
1 

x — — 

(771S2 - 4t/2M|. + 772m 2 + 4?7|M|, - m 2 t ) 

x {8m 6 Ms*(A; • p) 2 (p • p')(t ■ p') - 8m b M% c (m b + 2Mg, - m c )(k ■ p) 2 {t ■ p') 

+ 4M s *m c m 2 (ml + 4m b M s * + 4m c M s * - m 2 c )(k ■ p) 2 

- 8m b M B *(k -p)(k- p')(p ■ p')(t ■ p) - 32 (A; • p)(k ■ p')(t ■ p')m b M%* 
+ 8MB*m 2 (m b + m c )(k ■ p)(k ■ p')(p ■ p') 

+ 8m b M\ t (m b + 2M B . - m c ) (k ■ p){k ■ p') (t ■ p) 

- 16M|*m 2 (m 2 + m b M § , + m c M B * - m 2 c )(k • p)(k ■ p') 

- 8m b M§* (k -p){k- 1) (p • p') 2 + 32m 6 m c M§. (m b + M § * )(k ■ p)(k ■ t) 

- 8m 6 M§. (m b + 6M B , + m c ) (k ■ p)(k ■ t) (p • p') 
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Y?Re{M!M*) 



+ AM B . m 2 w {m b + m c )(k ■ p) (p • p')(t ■ p') 

+ 8M S * m c m 2 w (m 2 b + Am b M n , + 4Mj§ 

- 8M%,m 2 w (?>ml + hm b M%* + 2m b m c + 5M§*m c — m 2 c ){k • p)(t • p') 
+ ?>2m b M% t {k-p') 2 {t-p) - 32Ml i m 2 t (m b + m c )(k ■ p') 2 
+ 64m 6 M|. (M§. + m c ) (A; • p') (A; • t) 
+ 12M B *m 2 w (m b + m c )(fc • p')(P ' • p) 

+ 8M% t m 2 w {-m 2 b + m 6 Mg« + 2m b m c + m c Mg« + 3m 2 )(A: • p')(t ■ p) 

- UM%m 2 w {m b + m c )(k ■ p')(t ■ p') - 4M^m 2 w (m b + m c )(k ■ t)(p ■ p') 2 

+ 8M|.m* (mg + m b M B * c + m c Mg* - • *)(p ' p') 

- &M% t m c m 2 w (m 2 b + 4m 6 Ms* - 2m 6 m c + AM^m c - Sm 2 c )(k ■ t) 
+ \^m b M B ,m 2 w {p -p) 2 {t ■ p) - 8M B ,m 2 t m 2 w {m b + m c )(p ■ p') 2 
+ 16m b Ml ; m 2 w (m c + 2M B ,)(p ■ p')(t ■ p) 

- 48m 6 M|*ra 2 (p ■ p'){t ■ p) - 4Mf ; m 2 m 2 (-m 2 + 8m b M s * 
+ 8m c M s * + m 2 c )(p ■ p') - 32m b M^m 2 w (2m b + 3Mg*)(f ■ p') 

- 8m b m c M\ t m 2 w {m b — 8Mg* + m c )(t ■ p) 
+ lQm c Mg t m 2 rn 2 v (2rnl — m b Mg* + 2m b m c — m c Afg,)} 

In the above 

k ■ p = p ■ k = ^(s 2 - 4M| ; - m 2 w ), 
p-t~t-p = p-p' + p-k + AM L , 
k • p' — k • t — k ■ p — m 2 w , 

APPENDIX B 

In this appendix we give the expressions of squared matrix for t — > W + + b + T. The 
vertex coefficients are also neglected. 



p-p' = -{si-m 2 c -AMl t ) 

k-t = t-k = --(si - m 2 - m 2 ) 

* -P' = ~\( s 2 ~ m 2 - m 2 c ) 
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(s 2 + m 2 w - 2(m| + mj)) 2 {s 1 - 1m\ + m 2 b ) 2 



x {— • P ) 2 b • -p)-^r(k- P ) 2 (p - P ') + ^(k- P ) 2 (t ■ p) 

+ —( m b+ m t+^ m l)( k -p) 2 (t-p') 2 m t m b( k -P) 2 

m w m w 

+ 2m t + 2 <)( k -p)( k - p')(t ■ p) + 5^ • p)(* • k )(p ■ p 1 ) 

+ -^(-m 2 b + m 2 + 4m 2 ) (A; ■ p)(t ■ k) + 8(k ■ p)(p ■ p')(t ■ p) 

- 40m 2 (A; -p)(p- p') + 24m 2 (A; ■ p)(t ■ p) + 88m 2 b (k ■ p)(t ■ p') 

— 168m 2 m 2 (k ■ p) ^(m^ + m 2 — 2m 2 w ){k ■ p')(t ■ k) 

m w 

+ 24ml ( k ■ P')(t ■ P) - 32m 2 b m 2 t (k ■ p') + 24m 2 (t ■ k)(p ■ p') 

+ 96m^(t • k) + 4(2m 2 + 2m 2 - m 2 w )(p ■ p')(t ■ p) - Q4m 2 b m 2 (p ■ p') 

+ 16m 2 (2m 2 + 2m 2 - m 2 w )(t ■ p) + 8m 2 (7m 2 + 7m 2 t - 5m 2 )(t • p') - 17Qm b m 2 } 



x {8m 2 (p ■ p')(t ■ p) - 16m^(t • p') - 8m b (t • p) 
+ 8m 2 b (p-p')(t-p') + 2(p-p') 2 (t-p') 

- 3 -^(k.t)(k.p')- l -^(k.t)(k.p) 

m 2 w m 2 w 

+ -^{k-t){k.p'){p.p') 

m w 

1 fim 2 4 
+ —^-{k -t)(k ■ p)(p ■ p) H ^{k -t)(k ■ p')(p ■ p') 2 } 



(s 1 - ml) 3 (s 2 + m 2 w - 2{m 2 b + m 2 )) 



x {—^{k ■ P) 2 (P ■ P')(t ■p')-^4(k- pfit ■ p') + *^(* ■ p) 2 
m 2 w m 2 w m 2 w 

2 Am 2 

- —(k-p)(k-p')(p-p')(t-p) H ±(k ■ p)(k ■ p')(p ■ p') 

m 2 w m 2 w 
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m w m w 

- ^^(k .p)(k. p') -^(k.p)(k-t)(p. p'f 

ml m 2 w 

16ml,. Wl w 16mf., w , s 

JL(k -p)(k- t)(p ■ p') + —±(k -p)(k- t) 

ml ml 

+ 2(k ■ p){p ■ p')(t ■ p') + 16m 2 b (k ■ p)(t ■ p) - 2Sm 2 b (k ■ p)(t ■ p') 

+ s -4(k . m • p) - l -^ik . pr + ^4{k . P ')(k . t ) 

ml ml ml 

+ 6(k ■ p'){p ■ p')(t ■ p) + \2m\{k ■ p')(t ■ p) - 32m 2 b (k ■ p')(t ■ p') 

- 2(k ■ t)(p ■ p'f + Am 2 b {k ■ t)(p ■ p') - 8m 4 b {k ■ t) 

+ 4(p-p'f(t -p) -Am 2 (p-p'f + 12m 2 b (p ■ p')(t ■ p) 
12m\{p ■ p') (t ■ p') - 16m 2 m 2 (p ■ p') + 12m 4 b (t ■ p) 

- 40m$(t ■ p') + 8m 4 b m 2 } 



Here 



k-p = p-k=-(s 2 - Am 2 b - m 2 J, 
p ■ t = t • p — p • p' + p • k + 4ml, 
k-p' = k- t — k-p — m 2 w , 



p-p' = ^(si - 5m fe) 

k ■ t = t ■ k — —-(si — m 2 t — ml) 

t-p = --{s 2 -m t -m b ) 
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The Feynman diagrams for t 
Fig. 2 The branching ratio R 
Fig. 3 The branching ratio R 
Fig. 4 The branching ratio R 
Fig. 5 The branching ratio R 



Figure Captions 

B*W + c and t — > TW + b process at leading order in a, 
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